EXISTENCE OF STEADY SUBSONIC EULER FLOWS THROUGH 
INFINITELY LONG PERIODIC NOZZLES 



CHAO CHEN AND CHUNJING XIE 

Abstract. In this paper, we study the global existence of steady subsonic Euler flows 
through infinitely long nozzles which are periodic in x\ direction with the period L. It is 
shown that when the variation of Bernoulli function at some given section is small and 
mass flux is in a suitable regime, there exists a unique global subsonic flow in the nozzle. 
Furthermore, the flow is also periodic in x\ direction with the period L. If, in particular, 
the Bernoulli function is a constant, we also get the existence of subsonic-sonic flows 
when the mass flux takes the critical value. 



1. Introduction and Main Results 

The study on subsonic and transonic flows in nozzles has grown enormously in recent 
years. Subsonic and subsonic-sonic potential flows in infinitely long nozzles were studied 
in [201 ED HHJ EE] • For full compressible Euler equations, Xie and Xin in [22] showed the 
existence of global subsonic flow in an infinitely long nozzle which tends to be flat at far 
field. The key point in |22j is to transform the system of Euler equations into a second 
order equation of stream function. The careful energy estimates give far field behavior and 
uniqueness of flows. The idea in [22] was generalized to subsonic Euler flows in axially 
symmetric nozzles [14] . Subsonic and subsonic-sonic potential flows past a body were 
studied in [31 HI [7] and references therein. Subsonic Euler flows with nonzero vorticity in 
half space was investigated in JTU] . Subsonic flows were also studied as a part of stability 
of transonic shock problem, see [8j |9j O El [TlJ [121 El HI 1231 [26] and references therein, 
where subsonic flows and nozzles are small perturbations of some given background flows 
and nozzles with simple geometries, respectively. 

In this paper, we study the existence of global steady subsonic Euler flows through 
periodic nozzles. Consider 2-D steady isentropic Euler equations 

{pu) Xl + (f>v)x 2 = 0, (1) 
(pu 2 ) Xl + {puv) X2 + p Xl = 0, (2) 

(puv) Xl + {pv 2 ) X2 + p X2 = 0, (3) 
l 



where p, (u,v), and p = p(p) denote the density, velocity, and pressure, respectively In 
general, it is assumed that p'(p) > for p > and p"(p) > 0, where c(p) = yV(p) is called 
the sound speed. The most important examples include polytropic gases and isothermal 
gases. For polytropic gases, p = Ap 1 where A is a constant and 7 is the adiabatic constant 
with 7 > 1; and for isothermal gases, p = c 2 p with constant sound speed c [13J. 
We consider flows through an infinitely long periodic nozzle given by 

tt = {(x 1 ,x 2 )\fi(x 1 ) <x 2 < f 2 (x 1 ), -00 <x 1 < 00}, 

where /j (i = 1, 2) is L— periodic, i.e., fi(x\ + L) = fi{x\) for x\ G R. Suppose that there 
exist a G (0, 1) and C > such that 

\\M\c^ m < c and inf - /i(^i)) > 0. (4) 

si£[0,i] 

Therefore, the domain Q satisfies the uniform exterior sphere condition with some uniform 
radius r > 0. Without loss of generality, we assume that /i(0) = and /2(C)) = 1. 

Suppose that the nozzle walls are impermeable so that the flow satisfies the no flow 
boundary condition 

{u, v) ■ v = on d£l, (5) 
where v is the unit outward normal to the nozzle wall. It follows from ([T|) and (jSJ) that 

{pu, pv) ■ ndl = m (6) 



holds for some constant m, which is called the mass flux, where / is any curve transversal 
to the x\— direction, and n is the normal of I in the positive Xi-axis direction. 

Using the continuity equation (pQ), when the flow is away from the vacuum, the momen- 
tum equations (j5J) and (J3J) are equivalent to 

uu Xl + vu X2 + h(p) Xl = 0, (7) 
uv Xl + vv X2 + h(p) X2 = 0, (8) 

where h(p) is the enthalpy of the flow satisfying h'(p) = p'(p)/p and can be determined up 
to a constant. In this paper, for example, we always choose h(0) = for polytropic gases 
and h(l) = for isothermal gases. After determining this integral constant, we denote 
H = M p>0 h(p). 

It follows from © and (jSJ) that 

(u,v)-V{h{p)+ l -{u 2 + v 2 )) = Q. (9) 
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This implies that 



u 2 +v 2 
2 



hh(p), which is called Bernoulli's function, is a constant along each 



streamline. For Euler flows in the nozzle, we assume that at X\ = 0, Bernoulli function is 
given, i.e., 



where -80(^2) is a function defined on [0, 1]. 

When the Bernoulli function B is a constant, Proposition [3] shows that the flow is 
irrotational. The existence of periodic potential flows with small mass flux in periodic 
nozzles was obtained in [19J. In this paper, we first study the subsonic and subsonic-sonic 
periodic potential flows with relatively large mass flux. 

Theorem 1. If B (x 2 ) = B > H , then 

1. there exists an rh > 0, such that for any m e (0, m) there exists a unique subsonic 
periodic flow (p, u, v) which satisfies inf^ u > 0; 

2. the maximum of Mach numbers of the flows increases as m increases and goes to 
one as m — >■ rh, i.e., the flows approach sonic; 

3. there exist a sequence m n — > rh such that the associated potential flows (p n ,u n ,v n ) 
converge to (p, u, v) almost everywhere, which satisfies 



and the boundary condition (J5]) in the sense of divergence measure field, where p 
is determined by u and v via Bernoulli law. 

When Bq is not a constant, we have the following results on subsonic Euler flows in 
periodic nozzles. 

Theorem 2. Let the nozzle satisfy ^ and B in [TO}) satisfy 




(10) 



V x (u,v) = 0, 
div(pu, pv) = 0, 



(11) 



B' o (0) >0, B' (1)<0. 



(12) 



For any m e (0,rh), there exists eo > such that if Bq{x2) satisfies 



Bo — 5||ci,i([ ,i]) = e — e o> 



(13) 



where B is the constant in Theorem^ then 
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1. (Existence) there exists a periodic flow, i.e., 

p(xi + L, x 2 ) = p(xi, x 2 ), u(xi + L, x 2 ) = u(xi, x 2 ), v(x 1 + L, x 2 ) = v(x 1 , x 2 ), 

which satisfies the original Euler equations the boundary condition 

mass flux condition (TJ|) ; and the condition l[W\) ; 

2. (Subsonic flows and positivity of horizontal velocity) the flow is globally uniformly 
subsonic and has positive horizontal velocity. More precisely, there are two positive 
constants c\ and c 2 such that 

sup(u 2 + v 2 — c 2 (p)) < — c\ and u > c 2 ; (14) 
n 

3. (Regularity) the flow satisfies 

IIpIIc 1 ."^), ||w||c 1 ^(n)7 IM|ci.«(n) < c 

for some constant C > 0; 

4. (Uniqueness) the flow is unique in the class of the periodic flows satisfying (|T4l) . 

A remark about Theorem [2] is as follows. 

Remark 1. Using the analysis in this paper, it is easy to show that there exists a subsonic 
full compressible Euler flow in the nozzle, if the entropy is also prescribed at x\ = 0. 

The rest of the paper is organized as follows: In Section [2, we introduce the stream 
function formulation for the Euler equations and give the proof of Theorem [TJ In Section 
[31 a boundary value problem for stream function is analyzed. This is divided into two 
steps. The existence of solutions for the associated problem is studied in Section 13.11 The 
uniqueness, periodicity, and positivity of horizontal velocity of the flows are proved in 
Section 13.21 In Section HI we use the fixed point theorem to show the existence of Euler 
flows; the uniqueness of these flows are obtained by the energy method. 

2. Stream Function Formulation of the Euler Flows and Proof of 

Theorem Q] 

2.1. Bernoulli's Law and A New Formulation of the Euler Equations. We recall 
that the steady Euler system for subsonic flows is a hyperbolic-elliptic coupled system 
[22] . Therefore, one has to resolve the hyperbolic mode. 

To overcome the difficulties mentioned above, we introduce the stream function for 
2-D steady compressible Euler flows so that the Bernoulli function can be reduced to a 



single-valued function of steam function. This gives an equivalent formulation for Euler 
flows in terms of stream function. 

Proposition 3. For a smooth flow away from vacuum with no stagnation point, i.e., 
u 2 + v 2 > 0, the Euler system ([I])-© is equivalent to the system of equations (pQ) ; (191) , 
and 

vd xl B — ud xn B 
u 2 + v 2 

where B = \{u 2 + v 2 ) + h(p) and uj = d Xl v — d X2 u are Bernoulli function and vorticity, 
respectively. 



U = > ( 15 ) 



Proof: Let us first show that ([I])-© implies (fT5j) . Differentiating the Bernoulli function 
with respect to x\ and x 2 , respectively, gives 

d\B = udiu + vdiv + dih(p), d 2 B = ud 2 u + vd 2 v + d 2 h(p). (16) 

This, together with (J7j)-([H]), yields 

d\B = v(d 1 v — d 2 u) = vtu, d 2 B = —u(div — d 2 u) = —uuo. (17) 

Therefore, the equation (|T5|) holds provided u 2 + v 2 > 0. 

Conversely, it follows from straightforward computation that ([9]) and ( TT5l) imply ( iT7j ). 
Substituting (E]) into (EE]) gives ©and (El). Using ©, one has ©and 

This finishes the proof of the proposition. □ 

The continuity equation ([!)) implies that there exists a stream function ip such that 

d Xl ip = -pv, d X2 tp = pu. 
Hence, for the flows away from the vacuum, (J9j) is equivalent to 

■ VB = 0, 

where V -1 = (— d X2 , d Xl ). This yields that B and i[) are functionally dependent. Therefore, 
one may regard B as a function of ip. We denote this function by B = B(ip). It follows 
from the no flow boundary condition ©, that the nozzle walls are streamlines, so ip is 
constant on each nozzle wall. Taking into account, one may assume that 

ip = on 5*i, and ip = m on 5*2, (18) 

where S t = {(x 1; / 4 (xi))|xi 6 R} (i = 1, 2). 
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As in [22], for any given s > H , there exist g = g(s), g = g(s) and T = T(s) such that 

h(g(s)) = s, h(g(s)) + ^ = s, and p'(g(s)) = T 2 (s), 

where g(s), g(s), and T(s) are the maximum density, the critical density, and the critical 
speed, respectively for the states with given Bernoulli constant s. Set 

E(s) = Q( S )y/2( S -h(g(s))). (19) 

Then the straightforward calculations show that 

dp dg , dE 

> 0, > 0, and — > 0. 
as ds ds 

Obviously, g(s) < g(s), if s > H = mi p>0 h(p). Furthermore, there exists a 5 > such 

that 

0{si)>g{s2) for any s u s 2 e {B - 5,B + 5). (20) 
For a fixed s, if p and Ai satisfies 

Hp) + ^ = s, (21) 

then p is a two-valued function of Ai for Ai G (0, S 2 (s)) and the subsonic branch satisfies 
p > g(s), see [22]. When s varies, the subsonic branch will be denoted by 

p = H(M, s) for (M, s) G {(M, s)\M G (0, S 2 (s)), s > H }. (22) 

In view of (1211) . we have 

dH _ H 3 dH _ H 

~ds ~ H 2 c 2 -M > ' dM ~ 2(M - H 2 c 2 ) < ' 

2.2. Potential Flows and Proof of Theorem [0 If B (x) = B, we have B(ip) = B. 
It follows from (T22]) that p = H(\ Vip\ 2 , B). Furthermore, f|T5|) implies to = because 
Xi,X2) = -B in fi. Therefore, ?/> satisfies 

div ( „ =v ] = 0. (23) 

Let {M n } be a strictly increasing sequence satisfying lim n ^oo M n = Y?(B). We define 
H n (-,B) G C°°(M) satisfying 

H{M, B) if M<M„, 

H n (M,B)= { M E 2 (m (24) 

H(M n ,B) if M> n+ 1 j . 
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Combining Lemmas 2.1 and 3.1 in [20j and Lemma 1 in [19] (or Theorem 2 in [2]), there 
exists a unique periodic solution ip(-;t) of the problem 




ff„(|VVf,B), 
on Si, i/) = t on 

for t > 0. Define *JZ n {t) = sup^ |V"0 n (-; t)| 2 . Then that ^ n (t) is a continuous function of 
t follows from the same argument in Lemma 4.1 in [2U]. Let m n = sup{t|^# n (t) < M n } 
and rh = sup n m n . Then as m — > rh, the maximum of \Vip\ of solutions of the problem 
( |23l) and ( |T8l) tend to X(-B), i.e., the flows approach the sonic state. Moreover, Using the 
compensated compactness framework in [7] (Theorem 2.1) and [20] (Theorem 5.1), there 
exist a sequence {m n } f fh such that the associated u n = H ^® 2 ^™ 2 ^ , w n = — ^y^p §) ' 
and p n = 2 , -B) satisfy w n — )■ w, u n — > v, and p n — > p a.e., where (u, satisfies 

supf2 U pi\l\ — I; the system (Till , and the boundary condition ([5]) in the sense of divergence 
measure field. 

The proof of Theorem [1] finishes after we prove the following lemma on the properties 
of subsonic potential flows. 

Lemma 4. Ifm > 0, then the solution ip of the problem ( 12"3"|) and (|18[) satisfies inf^ dity > 
0. Furthermore, as m increases, maxjj |V"?/>| also increases. 

Proof: Note that ip satisfies < ip < m, therefore, ip achieves its minimum and maximum 
on the boundaries Si and 52 respectively. It follows from the Hopf lemma (Lemma 3.4 in 
[T6] ) that d X2 ip > on dVL. Therefore, the continuity of d X2 ip yields 

inf d X2 il> > 0. 

0<X!<L 

If ip satisfies f[2"3"l) . then d X2 ip satisfies 

T\\ n v, - (^ 2 c 2 -|V^| 2 )^ + ^M^ 
d l {a ij d J {d X2 ^)) = U, where a {j = (25) 

Here and later on, the repeated index means summation from 1 to 2. It follows from 
the strong maximum principle (Theorem 3.5 in [16]) that dity > on Cl. In fact, if the 
minimum of d X2 ^ is achieved at some point (x^xlj) G f2, by periodicity, we can always 
assume |a;*| < L. Then, it contradicts with the strong maximum principle for the equation 
(!25|) in the domain {(x 1: x 2 ) G Q : \x\\ < ^}. 
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It follows from Bernstein estimates (Theorem 15.1 in [16] ) and periodicity of the solu- 
tions that 

max|V?/>| < max|V , 0|. (26) 
n 9o 

Thus using comparison principle, Hopf Lemma and periodicity, we can show that the 
maximum of the flow speed increases as m increases (cf. Lemma 4.4 in [20]). □ 



2.3. Stream Function Formulation of General Euler Flows. For any m G (0, rh), 

there exists a unique periodic solution ^(xi,x 2 ) G C 2 ' a (Cl) for the problem (125]) and (ITS]) . 
Furthermore, there exist positive constants ctq and o\ such that 



< inf ip X2 = a <(J 1 = sup (27) 



Given W £ S defined by 



S = {WeC x ' p ([0,1]), [ W(s)ds = m, ||iy-^ 2 (0,-)||cM[ ,i] <^o/2}, 

Jo 

where (3 G (0,a), then W(s) > cr /2 for s G [0,1]. Therefore there exist a function 
y = k(i/j) such that 

1 f,= / (28) 



Differentiating ( 128]) with respect to ip yields 

W(k(V0) 

This shows that k(^) G C 2,/3 ([0, m]). Suppose that the Bernoulli function at X\ — is 
Bq(k(i(})) in terms of ip, since the Bernoulli function is a constant along each stream line, 
we get the Bernoulli function defined in the whole domain Q by B(x±, £2) = B(ip(xi,X2)) = 
B (K(if)(xi,x 2 ))). Define B(if>) = B (k(iJ))). Then B(ip) G C 2 ([0,m]). Combining (]IJ) and 
029]) gives 

||# - 5|| cl ,i([o, m ]) < Ce, B'(0) > and B'{m) < 0. (30) 

Therefore, p = H(\ Vip\ 2 , B(ip)) follows from (122]) for subsonic flows. The equation (IToT) 
becomes the following second order equation for the stream function if>, 

H mm*m) ) = m ^' K + m * ) - (31) 
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We first solve the equation f l3Tj) with the boundary condition (j!8j) . Second, we define 
a map from W to ip X2 {Q, •)• The fi xe d point of this map will induce the existence of the 
solutions of Euler equations. 

3. Analysis of the Boundary Value Problem for Stream Function 

3.1. Existence of Solutions. There are two main difficulties to solve the equation ( 13~TT) 
in Q. The first difficulty is that the equation (13TJ) becomes degenerate elliptic at sonic 
states. In addition, H is not well-defined for arbitrary ip and |V"0|; neither is B. The 
second difficulty is that this is a problem in an unbounded domain. Our basic strategy 
is that we extend the definition of B appropriately, truncate \Vip\ appeared in if in a 
suitable way, and use a sequence of problems in bounded domains to approximate the 
problem (I3T|) and (jlg|) . 



Set 



' B'{s), if < s < m, 



9{s) = 




It is obvious that g G C 0,1 (M) and 



^(^llc^OR 1 ) ^ l|£'( s )llc°'i([o,m]) < 2e/o- - 



Define 




Then, ||B'|| Cf o,i (R ij = ||5||cM(Ri) < ||B'||ci([o,m]) < Ce - Therefore, 



\B(ip)-B\ < Ce. 



Hence, there exists eo > such that if < e < eo, then B(ip) > H = inf s h(s). In view of 
fl3D|) . B also satisfies 



B\s) > for s < and B\s) < for s > m. 



(32) 



Let B = min^g^i] -Bo(^)- Choose 6q to be a fixed positive constant satisfying 



< 9 < min{E 2 (i?)/2, Y?(B) 




(33) 
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where a x is defined in ([27]). Let C G C°°(M) satisfy 



COO 



s, if s < -6*0/4, 

-0 o /8, if a > -9 /8, 



and define p = #(C(|V^| 2 - S 2 (i3(^)) + S 2 (S(^), <B(V0)- Instead of the equation (jgljl. 
we begin with investigating the equation, 

ft(^) + ft(^ fe )=pB'(V')- (34) 
The equation ( |34l) can also be written in the following non-divergence form 

A,(Vii>, ^dijip = V) (35) 

where 

C'p 

Ay(V^,^)) = p% + _ —-—- drfdjij), 

and 

'. 3 p 2 £ 2 -(C + S 2 ) + |V^| 2 , (C'-I)EE'IWP 



p 2 c 2 - (c + s 2 ) 

where the variables in ( and S are |V?/>| 2 — i3(^) and t3(ip), respectively. The direct 
calculation shows that the eigenvalues A and A of [Aj]2x2 satisfy < |A/A| < C for 
some constant C and thus the equation fl35|) is uniformly elliptic. However, J 7 involves a 
quadratic growth in so it is not easy to apply the classical elliptic theory directly. 

The strategy is to modify T by 

*W,*)=*W ^_^ +g) , (36) 
It is easy to see that J 7 = J 7 when ip satisfies |V"?/>| 2 — £ 2 (£>(?/>)) < —Oq/A. 

Proposition 5. There exists a solution ip(-;t) G C 2,Q! (f2) of the problem 

Aijdijip = f 

ijj = on Si, if) = t on S2. 
Furthermore, ift<m, then 

< ^ < m. (38) 
T/iere exzsfo mi > and e\ > suc/i t/iat if < t < mi and \\B'\\ < e x , then 

|V</f - £(W)) <-V2- (39) 



(37) 
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Proof: For the unbounded domain Q, we use a sequence of boundary value problem 
defined in bounded domains Qn to approximate the problem in Q, where Qn £ C 2,a 
satisfies Qn C C f^iv with ttj. := R < A;}. The construction of can be 

found in (20]. 

We first solve the boundary value problem 



A ij (Vil),il))d ij il)=J r (Vil),il)) in tt N , 



X2~fl(xi) 



-t on dQiv. 



(40) 



/a(:ci)-/i(a!i) 

Similar to the proof of Proposition 3 in [22], there exists a C 2,Q -solution ip^, such that 



\^n\<C [t + 



T_ 

A 



< C A/A, |/i|c 2 .«,m, 



A 



By Arzela-Ascoli theorem, one can select a subsequence of {4>n} (We still label it by 
{4>n}), such that 

ip N -> V> in C 2 ' ai (^) with KiO, ai < a. 



Obviously, t/> solves (]37|). 

Since £> satisfies (I3"2"j) . then ^(ViPn^n) > in the domain fl N D {^jv > m }- Thus 
according to system (j4"0~]) . 



Aij(Vi/)N,i*N)dijil)N > in n {^at > m}. 

Since 

by the maximum principle (Theorem 3.1 in [16]), one has 

ipN < sup?/^ < m in D {V'at > m}. 

Similarly, it is also true that 

ipN > inf ipN > in Qn H {^iv < 0}. 

c9f2jv 

Therefore, one has < ipN < m. So the limit ip satisfies ( 138|) . 

The Holder estimate for the gradients of elliptic equations of two variables implies that 



|W>7v(x) -W>iv(j/)| . ... . , , . , .. , 

Pli,^, = sup — — < C(A/A, |/i| 2 ) ( 1 + t + 



\X — y\P 
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Then, ip^ satisfies the estimate 

Wn\ 2 <v\ i + t- 



c v I t 



T_ 
J 



where C v is independent of N. Note that 

\F\ < Ce 

there exist rji, mi and e\ such that 

771 (1 + m\ + Cti) 2 < min{ 
If t < mi, e < 6i, then 



" 4^f } ' ^( mi + C ei ) 2 <m.n{^i,j}. 



|Wjv| 2 -S 2 (B(^)) < 



< 



Thus, the limit ip satisfies ( 139]) if ei is suitably small. 
This finishes the proof of the proposition. 



(41) 



(42) 



□ 



3.2. Uniqueness, Periodicity and Positivity of Horizontal Velocity. In this sub- 
section, we show that the uniformly subsonic solution obtained in Proposition [5] for the 
problem ( l3Tj) with the boundary condition 

■0 = on Si, and if) = t on S2. (43) 

is unique. 

Proposition 6. There exists a positive constant €2 < €1 such that if \\B'\\co,i = e < e 2 , 
then for < t < m, the solution if) of l[JJ\) and (I43p which satisfies 

|W| 2 - E(B(V>)) < -^o/4, < i) < m (44) 

must be unique. 

Proof: Let ipi,i = 1,2, both solve (EE)) and g3]). Set * = Vi - ^2- Then ^ satisfies the 
following elliptic boundary value problem, 

d^d^) + = + <M, in Q, 

* = 0, on dfl, 
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where 



7o (|W| 2 -F 2 c 2 )# ' ' " Jo # 2 c 2 -|W| 2 

7o |W| 2 -# 2 c 2 7 # 2 c 2 -|V^| 2 



with 5 = Jj/(H), H = H(\W\ 2 ,B(^)), and # = + (1 - 6)ip 2 . Since ||£'||co,i < Ce, 
we have 

I foil + hi + \d\ < Ce. 
Choose a smooth cut-off function r\ £ C^°(R) satisfying 

„ = „(.)=(<>■ I'l^^* (46) 
[1, \s\ <NL. 

Multiplying r] 2 (xi)^ on both sides of the equation (T45|) and integrating by parts yields 
that 



/ / \V^\ 2 dx 1 dx 2 <C r] 2 (a ij d i ^d j ^)dx 1 dx 2 

J J{\xi\<NL}nU J Jn 

< C JJ (|W*I(I*I + IV* |) + e^ 2 |^|(|^| + |V*|))tfaida; 2 



< C // (|^| 2 + |W| 2 )c£Mx 2 

+ Ce [[ (\V\ 2 + \W\ 2 )dx 1 dx 2 . 

J J{\xi\<NL}nQ 

Since \1/ = on <9f2, by Poincare's inequality, we have 

\^\ 2 dx 1 dx 2 < C II \V^\ 2 d Xl dx 2 

{fcL<|xi|<(fc+l)L}nn J </{fcL<|zi|<(fc+l)L}nn 

for any fceZ. When e < e 2 for some e 2 > 0, we have 

\W\ 2 d Xl dx 2 <C II \W\ 2 d Xl dx 2 ., (47) 

{\x!\<NL}nQ J J{NL<\xi\<(N+l)L}nn 

Noting that |W| G it yields that W G L 2 (fi). Thus 

lim / / |W| 2 Gbidx 2 = 0. 

N ^°° J J {NL<\x!\<{N+l)L}nn 

Therefore, the estimate (T4T1) implies that jj n |V\I/| 2 = 0. It follows from $ = on dQ 
that = in Q. Therefore, the solution of fl3Tj) and fj43l) is unique. □ 
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One can check easily that if ip(xi,X2) solves the boundary value problem f l3T|) and (fT8l) . 
so does ip{x\ + L,X2). Then the uniqueness implies the following corollary. 

Corollary 7. For any 6q > 0, there exists a positive constant €2 < ex such that if 
\\B'\\co,i = e < €2, then for < t < m, the solution ip of (3l\) and ( 1431) satisfying 
(|44|) must be periodic with respect to X\ with period L, i.e., 

4>(xi, x 2 ) = ip{xi + L, x 2 ), V(xi, x 2 ) G fi. 

Now we show that "^(O, •) G S. 

Proposition 8. For any 6 > 0, there exists a positive constant e 3 < €2 such that if 
||i3'|| c o,i = e < e 3 , then for < t < m, the solution ip of l[3l\) and fj43|) satisfying 
then \\Vip - V^llc 1 ." < Ce. 

Proof: Set = tp — -0. Subtracting ff23l) from fl3Tl) gives that \l/ satisfies 

^(a^* + 6<(B(V) - &)) = #(| W| 2 , Bty))^), in SI, 
* = on a 



(48 



where 



' #(# 2 c 2 - I W| 2 ) Jo H 2 c 2 - | W| 2 



with c = yp'(H), H = H(\\7^\ 2 ,6B(ip) + (1 - 9)B), * = 9ip + (1 - 0)0. Since both 
and are periodic, \I/ is also periodic with period L. Multiplying the equation in (1481) 
with \1/ and integrating the resulting equation on Q R {— 2L < X\ < 2L}, and integration 
by parts yield 

II V 1 I / || L 2( nn {_2L<xi<2L}) < Ce. 

Applying Moser's iteration (Theorems 8.17 and 8.25 in [IS]), we have 

ll*llL°°(on{-| <j:i<f}) - Ce. 

Using the estimate for elliptic equation of two variables (Theorem 12.4 and global esti- 
mates on page 304 in [16]), we have 

ll^llci.^nn-tCKaji^L}) < Ce. 

Using Schauder estimate, we have H^H^am^ < Ce. Choosing e 3 > sufficiently small 
shows that || V\I / ||c ,1 > t *(Q 1 ) < 0"o/2 provided < e < e 3 . This finishes the proof of the 



proposition. 
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□ 



4. Existence and Uniqueness of the Euler Flows 

In this section, We prove the existence and uniqueness of subsonic solutions for Euler 
equations. 

Proposition 9. There exists a positive constant e§ such that if \\B' Q \\ci < then the 
system ([I])-© under the conditions (TJP, mass flux condition fifij), and the condition (T77J) 
has a subsonic solution with positive horizontal velocity. 

Proof: It follows from Proposition that there exists a solution ip(-',t) for the problem 
d37D. Set 

S{t) = {ip(-,t)\ip(-;t) solves (J37D} 

and 

A(t)= inf 8up{|VV(-;*)r-S 2 W(-;t)))} 

V>es(t) q 

Let fh = sup s {s G (0, m)|A(s) < — /A}. Then A(t) is continuous for t G [0, m](cf. 
Proposition 6 in [22]). We claim that fh > m. Indeed, it follows from Proposition [5] that 
A(mi) < —6 /2. If fh < m, then 

|V^(-;m)| 2 -S 2 (i3(^(-;m))) 

= | V^(-; m)| 2 - |V$(-; m)| 2 + |V^(-; m)| 2 - £ 2 (5) + T?(B) - S 2 (i3(^(-; m))) 

<Ce -6 + Ce, 

where we use the property that the maximum of flow speed increases as the mass flux 
increases in Lemma HI Thus there exists a positive €4 < €3 such that |V^(-;m)| 2 — 
S 2 (i3(^(-; 'hi))) < —#o/2. This contradicts with A(m) = —9q/A which follows from the 
continuity of A(s). Thus fh > m. Thus it implies that for t = m the problem (13 7p has a 
unique solution satisfying 

|W| 2 -x(£(^)) <-b q /a. 

Using Proposition [HI there exists a positive constant €4 < €3 such that \Vip — V^| < cr /2. 
In particular, \d X2 ip(0,-) — d X2 ^(0, -)\ < cr /2. Therefore d X2 i/j(0, •) G 5. Hence, we can 
define a map T : S — > S by 

T(LU) = ^ 2 (0,-). 
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By Proposition El ^^(O, •) — ^^(O, -)||ci.« < Ce. Thus TS is a compact subset of 
S. It is easy to see that T is a continuous map. Note that S is a closed convex set in 
C 1,/3 [0, 1]; the existence of fixed point of T on S follows from Schauder fixed point theorem 
(Theorem 11.1 in [To]). 

Let W be the fixed point of T and k(y) satisfies 

rKv) 

y= I W{s)d.S. 

Define B(if)) = B (k(if})). Then 

if) = on Si and if) = m on S2 

has a solution satisfies "^(O, £2) = W / (a;2). It is clear that p = H{\ V^| 2 , B{ip)), u 
d X2 ip/p, and u = —d X2 ip/p satisfy the Euler equations and the condition 







+ h(p) (0,x 2 ) = B{if}{0,x 2 )) = B {x 2 ). 



u 2 + V 2 



Thus we get the existence of solution for Euler system under the conditions (JSJ), mass flux 
condition fl6]), and the condition ( flOl) . Furthermore, it follows from Proposition |8] that 
horizontal velocity is positive. □ 

Now we show that the periodic Euler flow is also unique. 

Proposition 10. There exists a positive constant such that if \\B f Q \\cn < 6$, then 
the uniformly subsonic solution the system ([I])-© under the conditions |3P, mass flux 
condition fifjj), and the condition (flQJj is unique. 

Proof: Suppose that ipi (i = 1,2) are stream functions of two periodic solutions of 
the Euler equations with positive velocity. Let (i = 1,2) be the functions satisfying 
Ki(if)i(0, X2)) = X2 (i — 1,2). Set Bi(if>) = Bo(ni(if})). Then ipi satisfies the problem 

div ( Trn ^ V f„ — )= H(\Vij\ 2 ,B t (iPm(iP), in O, 



.^(|V^| 2 ,^W), 
if> = on Si, if) — m on 5*2. 



17 

Set — ipx — ip2- Then \1/ satisfies the problem 

diia.jdj^) + diipiV) = Cidi^f + dV + e£, in fi, 

(49) 

^ = on dtt, 



where 



(H 2 c 2 - 


W| 2 )c5 


ij + i)M'<)/\> 


H{H 2 c 2 






2 ) 



a l3 = / v ' *' ' L~Z ' M, h = - I * ,._ rfg 



1 

ifV _ | V ^|2 



Jo |VW|^ — H 2 c 2 

d = f u~2 2 H3 ^tJ 9B '^ + WW*)* ' e = / HdB, 
Jo H 2 c 2 - |Vvp Jo 

with * = 0^1 + (1 - 0)V>2, # = H{\ V§| 2 , + (1 - flWVfe), p = WO - Mfo) 
and£ = ^(^i)-i3 2 (^2). 

We first multiply ^ on both sides of the equation in (j49l) and integrate the resulting 
equation on In view of periodicity of the coefficients and \I/ in (j4"9~|) . integration by 
parts yields 



(50) 



Thus, 

j Qi \VV\ 2 dx < Cej Qi {\y\ 2 + \Vy\ 2 ) + Cj ni {\V\ 2 + \£\ 2 )dx. (51) 

Since ^ G M/ 1 ' 00 , the first term on the right hand side of (15T|) is uniformly bounded. 
Note that \l/ = on dQ D fii, Poincare inequality implies that 

\^\ 2 dx <C I I W| 2 cfe. 

Therefore, the first term on the right hand side of (15T|) can be absorbed by the left hand 
side. 

Now let us estimate the second term on the right hand of flBTj) . First, 



/ \V\ 2 dx= |B °«i(^i) - B o K 2 (ij 2 )\ 2 dx 



< \B o «i(V>i) - 5 o ki(^ 2 )| daj + / l^o o Kl (^ 2 ) - So ^2(^2)1^ 
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Using mean value theorem, we have 

h<Ce [ \^\ 2 dx <Ce I \V^\ 2 dx. 

JUx JO! 

Note that ipi(0, Kity?)) = ^2(0, ^2(^2)), we have 

d X2 ipi(0, s) - d X2 i/j 2 (0, s)ds = / d X2 ip 2 (0, s)ds. 
In view of the fact that d X2 ip 2 > , we have 

|«i(^)-« 2 (^)l <C||V*(0,.)||^[o,i] 

Thus, 

h < J Ql Ce|| W(0, -)\\i?\p,i]dxidx2 < Ce\\ WH^m). 

Similarly, we can show that f n \£\ 2 dxidx 2 < Ce|| V^Hloo^^. 
This implies that 

l|V^|| L 2 (ni) < Ce||V*|| L «. (ni ). 
By Nash-Moser's iteration, we have 

00 (Oi) 

Using the estimate for elliptic equation of two variables (Theorem 12.4 and global esti- 
mates on page 304 in [IS]), we have 

ll^llc^Cfii) < ^11^11^00(0^ = Cell^Hioo^) 

Therefore ^> = in Q. 

This finishes the proof of the proposition. □ 

Choosing eo = min{ei, e 2 , 63, 65, e^}, then Theorem [2] follows from Propositions [9] and 
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